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Abstract 

The cross sections for the production of hadrons in quasi-real photon- 
photon collisions in proton-proton and heavy ion interactions are compared 
with the corresponding cross sections for central diffraction and for photon- 
pomeron collisions. The signatures, heavy ions or protons with only slightly 
changed momenta together with two large rapidity gaps and a cluster of 
produced hadrons in the central region, are nearly identical in all three 
processes. Therefore, it will be rather difficult to distinguish the reactions 
experimentally. It is found, that central diffraction is the dominant process 
in collisions of protons, light and medium-heavy ions. The photon-pomeron 
and photon-photon processes have quite similar cross sections in collisions 
of heavy ions like lead. 
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1 Introduction 



Double-photon exchange as well as double-pomeron exchange and photon- 
pomeron reactions are characterized by large rapidity gaps separating the rem- 
nants of the colliding hadrons or heavy ions from the particles produced in the 
central rapidity region. This feature is likely to give a considerable reduction of 
the background to the studied reaction. 

Photon-photon collisions in proton-proton or heavy ion reactions have been 
discussed repeatedly ]T|, 0, [| £|, || as a very attractive reaction channel. Processes 
studied include Higgs production || |7|, [|, ||, |H| and production of SUSY particles 
Also the ALICE Collaboration discusses this option for the heavy ion 
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experiment at the Large Hadron Collider (LHC) of CERN | 12fl . 

Central diffraction (CD) and especially hard central diffraction is recognized 
as an interesting new reaction channel at future hadron and heavy ion colliders 
at least since the first model for diffractive hard scattering due to Ingelman and 
Schlein [|13| and the pioneering experiment of the UA8 Collaboration [14], 15]. For 
example, Higgs production in this channel was discussed in Refs. 
and the production of heavy flavors was studied in Refs. 



16, 17, 18, 19 



20, 21, 22, 23 . The 



resonance and glueball production in photon-photon collisions and in central 
diffraction was compared in Refs. [^, . 

Here we will discuss both processes and in addition photon-pomeron collisions 
for the same reactions in the framework of the two-component Dual Parton Model 
(DPM) p6| , |27| , |28 | and in particular in the version as implemented in the DPM 
event generator Phojet ]29|, |30| . We restrict the paper to the calculation of 



cross sections for hadron production in the three processes, in a later step we 
might extend the calculation to cross sections for heavy flavor or new particle 
production. 

Hadronic photon-photon collisions within this framework were already studied 
and compared to data in a recent paper by Engel and Ranft |3D]] mainly for 
photon-photon collisions at electron-positron colliders. Here we report about the 
extension of the Phojet model to photon fluxes as expected at proton-proton 
and heavy ion colliders. 

Hard diffraction within the two-component DPM was first studied and com- 
pared to UA8 data and data from HERA by Engel, Ranft, and Roesler plf. In 



the present paper we will calculate central diffraction cross sections as well as 
photon-pomeron cross sections for proton-proton and heavy ion collisions. Using 
these cross sections we compare basic features of the three reactions, double- 
pomeron, double-photon and photon-pomeron interactions at proton-proton and 
heavy ion colliders, which have very similar experimental signatures. A second 
more detailed study of particle production in double-pomeron and double-photon 
scattering, also using the Phojet event generator will be published soon |32j| . 

In Section 2 we describe the models used. Section 3 gives a detailed compari- 
son of the cross sections and hadron distributions expected in the three channels. 
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A summary is given in Section 4 and in the Appendix we collect all the details 
about photon flux calculations in heavy ion reactions as used in this paper. 



2 The Models 

2.1 The event generator P HO JET 

The realization of the DPM with a soft and a hard component in Phojet is 
similar to the event generator Dtujet-93 J27I, [£J simulating pp and p-p collisions 
up to very high energies. Phojet is applicable to collisions of stable hadrons as 
well as of photons. Here we give only a short summary on the physics of PHOJET, 
for more detailed descriptions we refer to J29|, [30], [34], ^ . 



The interactions of hadrons and the hadronic fluctuations of the photon are 
described within the Dual Parton Model in terms of reggeon and pomeron ex- 
changes. The physical photon state is treated as a superposition of a "bare pho- 
ton" and virtual hadronic states having the same quantum numbers as the pho- 
ton. For soft processes, photon-hadron duality is used. The energy- dependences 
of the reggeon and pomeron amplitudes are assumed to be the same for all 
hadronic processes. Therefore, data on hadron-hadron and photon-hadron cross 
sections can be used to determine the parameters necessary to describe soft 
photon-photon interactions. Inelastic interactions are subdivided into processes 
involving only soft processes and all the other processes with at least one large 
momentum transfer (hard processes) by applying a transverse momentum cutoff 
^cutoff £ partons. On Born-graph level, for example, the photon-photon cross 
section is built up by: (i) soft reggeon and pomeron exchange, (ii) hard double- 
resolved photon-photon interaction, (iii) hard single-resolved interactions, and 
(iv) hard direct interactions. The Parton Model calculations of the hard pro- 
cesses have been done using the leading order GRV parton distribution functions 
for the proton |3T| and the photon 



The amplitudes corresponding to the one-pomeron exchange between the 
hadronic fluctuations of the photon are unitarized applying a two-channel eikonal 
formalism similar to ^7|. In impact parameter representation, the eikonalized 



scattering amplitude for resolved photon interactions has the structure 

a tm (s,B)= % -(^] 2 (l-e-^ B A (1) 



with the eikonal function 

X (s, B) = X s(s, B) + X h(s, B) + X d(s, B) + X c(s, B). (2) 

Here, Xi{ s iB) denote the contributions from the different Born graphs: (S) soft 
part of the pomeron and reggeon, (H) hard part of the pomeron (D) triple- and 
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loop-pomeron graphs, (C) double-pomeron graphs. In PHOJET the last two terms 
are included in the unitarization whereas in Dtujet they are taken in lowest 
order. The amplitude for hadron-hadron interactions has the same structure as 
(HD but without the couplings ej f q q of the photon to the hadronic fluctuations. 

The probabilities to find a photon in one of the generic hadronic states, the 
reggeon and pomeron coupling constants, and the effective reggeon and pomeron 
intercepts cannot be determined by basic principles. These quantities are treated 
as free parameters and determined by cross section fits [29]. In Refs. [f29, |3(J 



the model predictions for the proton-proton, photon-proton, and photon-photon 
cross sections are shown and compared to data. 

The probabilities for the different final state configurations are calculated 
from the discontinuity of the scattering amplitude (optical theorem) which can 
be expressed as a sum of graphs with k c soft pomeron cuts, l c hard pomeron cuts, 
m c triple- or loop-pomeron cuts, and n c double-pomeron cuts by applying the 



Abramovski-Gribov-Kancheli cutting rules |38], In impact parameter space 
one gets for the inelastic cross sections 



a(k c ,l c ,m c ,n c ,s,B) 



(2 Xs ) kc (2 X h) u (2 XD ) mc (;2\ ' 



exp[-2 X (s,B)}. (3) 



k c \ U m c \ n c 

Since the triple-, loop-, and double-pomeron graphs are objects involving sev- 
eral pomerons, a further resummation is done [f27| , |35f to allow the probability 
interpretation of Eq. (|3|). 

For pomeron cuts involving a hard scattering, the complete parton kinematics 
and flavors/colors are sampled according to the Parton Model using a method 
similar to pEOfl , extended to direct processes. For pomeron cuts without hard 
large momentum transfer, the partonic interpretation of the Dual Parton Model 
is used: photons or mesons are split into a quark- ant iquark pair whereas baryons 
are approximated by a quark-diquark pair. The longitudinal momentum fractions 
of the soft partons are given by Regge asymptotics PTl O. The transverse 



momenta of the soft partons are sampled from an exponential distribution in 
order to get a smooth transition between the transverse momentum distributions 
of the soft constituents and the hard scattered partons. 

In diffraction dissociation or double-pomeron scattering, the parton config- 
urations are generated using the ideas of the two-component DPM applied to 
pomeron-hadron, photon-pomeron and pomeron-pomeron scattering processes 
(see |31| and references therein). For the parton densities in the pomeron we 



use the CKMT parametrization [[43 



Finally, the fragmentation of the sampled partonic final states is done by 
forming color neutral strings between the partons according to the color flow. In 
the limit of large number of colors in QCD, this leads to the two-chain configu- 
ration characterizing a cut pomeron and a one-chain system for a cut reggeon. 
The chains are fragmented using the Lund fragmentation code Jetset 7.3 pH . 
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2.2 Diffraction dissociation in the two-component Dual 
Parton Model 



2.2.1 Single diffraction dissociation 

In the model, low- and high-mass diffraction dissociation is distinguished. Low 
mass diffraction dissociation is described by assuming the scattering of a superpo- 
sition of resonances with the quantum numbers of the dissociating particle [[15 



For simplicity, only one generic resonance is considered in the model. In this 
case, the cross section of diffractive processes with low-mass final states (Mf, < 2 
GeV 2 /c 4 ) can be calculated using a two-channel eikonal formalism [|6|, j29fl . 
In the limit of large diffractive masses s ^> Mjj ^> sq and M5 ^> t with sq ~ 1 
GeV 2 , the data can be understood in terms of the triple-pomeron graph ||6|, ^7[| . 
Here, Mp denotes the diffractively produced mass and t is the squared four- 
momentum transfer. In Fig. [l], the triple-pomeron graph and the corresponding 
diffractive cut for diffractive dissociation of particle A and quasi-elastic deflection 
of particle B is shown. Assuming multiperipheral kinematics of the pomeron cut 
final states, the rapidity gap is approximately given by 77 gap ~ ln(s/M|,). It is 
convenient to characterize the final state using the Feynman xp of the elastically 
scattered particle B 

which leads to 

M 2 D » (1 - \x B \)s and r/ gap ^ln(- — ,) (5) 

\l — \xb\J 

In Born-graph approximation neglecting rescattering effects, the double- 
differential cross section for high-mass diffraction dissociation reads 



d 2 a T p 1 / \ 2 ( s \ 2Af ( s o 



dtdMl 16vr 



n (qaip) 9°ipipip 9 bp 




a F (0) 



xex P {(^ + W + 2^(0)ln(^jj tj . (6) 

with Ajp = «jp(0) — 1, ajp(O) beeing the pomeron intercept. The coupling con- 
stants are parametrized by 

9ip{t) = g&expfybiptj i — A,B (7) 



gpipp(ti,t2, fa) — S'jpppexp (^-bpppfti + t 2 + t 3 )j 



(8) 
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In the following we will restrict the diffract ively produced mass according to 
an experimentally motivated cut on the Feynman xp of the elastically scattered 
particle |^0| 

\x B \>c, M 2 DMn <M 2 D <(l-c)s, (9) 

with c = 0.9 . . . 0.97. The integration over t and Mjj can be performed analyti- 
cally 



0"TP 



x 



X 



167T 

exp < 



9ppp 9 bp 
2<(0) 

bAP + bjpjpjp I 



. s 



2a' F (0) 



Ei (Ajp f ^ + ^ +ln 



2^(0) 




(10) 



where Ei denotes the second exponential integral function. 

Elastic and inelastic rescattering effects decrease the cross section as given 
by Eq. ([]) considerably |SD[ . This suppression is estimated using the eikonal 
model as described in the previous section. In impact parameter representation, 
the experimentally observable cross section of diffraction dissociation follows from 
(see Eq. (§) 



a d (s,B) = a(k c 



0,/ c 



0, m c — 1, n c — 0; s, B) . 



(11) 



The free parameters of the model (coupling constants, pomeron intercept, 
slope parameters) are determined by a global fit to data on total, elastic and 
diffractive cross sections as well as data on elastic slope of pp, pp and 7p interac- 
tions [p9|. Assuming that soft hadronic interactions in hadron-hadron, photon- 



hadron, and photon-photon interactions can be described by the exchange of one 
universal object, the pomeron, all these data can be combined to increase the 
predictive power of the model. 

Fig. [I] suggests the interpretation of photon diffraction in terms of photon- 
pomeron scattering. However, since the pomeron is only a theoretically intro- 
duced object to describe the some features of hadronic high-energy scattering, it 
is not possible to consider photon-pomeron scattering without the corresponding 
hadron where the pomeron couples to. 

Finally, it should be mentioned that the model results can be compared with 
the new measurements on the cross section of photon single diffraction dissoci- 
ation at HERA. This is important since the cross section on photon diffraction 
dissociation enters directly the predictions on the photon-pomeron cross sections. 



For 



200 GeV, the model predicts a cross section of 19/ib to be compared 
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with 23 ± ll//b Furthermore, in Ref. | ^T| , some model results are shown to- 
gether with data on photon diffraction obtaind by the HI Collaboration, finding 
reasonable agreement. 



2.2.2 Central diffraction 

A direct consequence of the interpretation of diffraction as pomeron-particle scat- 
tering is the prediction of the existence of pomeron-pomeron scattering (double- 
pomeron scattering). Double-pomeron cross sections in hadron-hadron collisions 
were first calculated on the basis of single diffractive measurements and Regge 
theory by Chew and Chew [|52] and by Kaidalov and Ter-Martirosyan Using 



different models, these calculations have been continuously extended to partial 
cross section estimations (see for example Refs. |55[| and references therein). 

The calculation of double-pomeron cross sections is subject to considerable 
uncertainty. There are several reasons for this: 

(i) The cross section is proportional to the square of the triple-pomeron coupling 
constant gjpipip, which is not very well known. 

(ii) The double-pomeron cross section differs considerably, if this term is included 
into an unitarization procedure for all hadronic cross sections or not. For exam- 
ple, it was discussed recently pBfl , that the effects of shadowing decrease the 
double-pomeron cross sections at LHC energies by a factor of 5. The cross sec- 
tion estimates as given for example in |53, 20] do not include any unitarization 



whereas in our approach the double-pomeron cross sections are included in the 
unitarization. 

(iii) The energy dependence of the double-pomeron cross section might differ 
considerably when using either a critical or a supercritical pomeron intercept 



ctjp(0). In |53], |2(| a critical pomeron intercept ajp(0) = 1 was used. Here we use, 
consistent with high-energy cross section measurements, a supercritical intercept 
a F (0) > 1. 

(iv) On a more practical level these cross sections depend on the cuts applied to 
the centrally produced cluster of particles and on the rapidity gaps demanded by 
the experimental triggers. 

(v) In CD, one or both of the incoming hadrons can be excited to a resonance 
N*. In the following, for the comparison with the two-photon cross sections, we 
exclude from the calculation all cross section contributions involving resonances. 
Including these contributions would about double the cross sections given below. 
These contributions would also change the rapidity gap demanded in experimen- 
tal trigger conditions. 

In the approximation of multiperipheral kinematics, the mass M c d of the cen- 
trally produced diffractive system is M^ d = (s\S2)/s (see Fig. ^|). The rapidity 
gaps between the central system and the elastically scattered particles can be 
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approximated by 

ffeapl ~ 111 [j-J %a P 2 « 111 . (12) 

Denoting the Feynman xp of the elastically scattered particles A and I? by xa 
and ib one gets PU| 

= 1 - - z B = - (l - -) (13) 



.S 



and 



M c 2 d «(l-|xAi)(l-|x B |) S 

W - In (j^) %a p2 - In ( T -^- T ) . (14) 

Within the framework of Gribov's Reggeon calculus, the amplitude of the 
double-pomeron graph can be calculated (for more details see j3*2"|). From 
this one gets the cross section as function of tA, ts, «i, and s 2 (tA and t B denote 
the squared momentum transfer of particle A and B). After integration over tA 
and ts the differential cross section reads 

d^dM[ d ~ 256^^ (MJ {M[J MfJ[ 



x 



(9°ajp) 



2 



si 



&af + bjpjpjp + 2a']p(0) In 
{9° B]P f 



b>BTP + ^iP/PJP + 2a'jp{Q) In 



(15) 



with 



2 / M 2 \ Ajp 

^(M C 2 J = (g° m ) (-^) . (16) 



Applying the previously discussed cut on the Feynman of the elastically 
scattered hadrons 

x A > c, \x B \ > c, 

M 2 cdMn < M 2 cd < (1 - c) 2 S (17) 
the integration over Si can be performed 

cd < si < (1 - c)s (18) 



c 
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1 cd 

1 / \ 2A JP 



<(0) ^M C 2 J 
xln / &AiP + &iP7PiP + 2a^(0) ln((l - c)s/M 2 d 
X n V b B1P + bjpjpjp + 2^(0) ln(l/(l - c)) 

x (&AiP + &BiP + 26jpjpjp + 2a^(0) ln(s/M 2 d ) 



(19) 



In order to calculate the cross section for CD, the Born graph cross section (19) for 
double-pomeron scattering is included in the eikonalization. In impact parameter 
representation, the CD cross section reads (see Eq. (Q)) 

a cd (s, B) = o(k c = 0,l c = 0, m c = 0,n c = 1; s, B) . (20) 

In Fig. | we compare as function of the energy the CD cross sections in proton- 
proton collisions, which we obtain from Phojet with the cross section obtained 
by Streng ||20|| . For both calculations the same three kinematical cuts are used: 
M cd >2GeV/c 2 and c =0.90, 0.95 and 0.97. In Phojet we use a supercritical 



pomeron with Ajp = 0.08 whereas Streng |20[ uses a critical Pomeron with Ajp 
= 0. Note that the double-pomeron cross section grows in Born approximation 
with s like ~ s 2Ap . This rapid increase is damped in Phojet by the unita- 
rization procedure. At high energies, contributions from multiple interactions 
become important. The demanded rapidity gaps are filled with hadrons due to 
inelastic rescattering and the cross section for CD gets strongly reduced. In con- 
trast, Streng calculates only the Born term cross section. Figure |^ illustrates the 
differences obtained using different methods. We stress, both methods use the 
measured single diffractive cross sections to extract the triple-pomeron coupling. 

Finally, it should be mentioned that in a very recent work j57| pomeron- 
pomeron cross sections have been studied in the perturbative BFKL-Bartels ap- 



proach (see [57| and references therein). The predictions on cross sections found 



in this work are of the same order as the ones obtained here. 



2.3 Diffractive cross sections in collisions involving nuclei 

2.3.1 Central diffraction cross sections in heavy ion collisions 

There are certain difficulties in resolving the CD cross section as a function of 
its central mass M c & in hadron-nucleus and nucleus-nucleus collisions. In fact, 
consider the simplest case of pd scattering with only two possible inelastic pN 
collisions. If in the first collision a central mass Mi is produced, in the second 
collision another central mass M 2 is produced and their rapidities overlap, then 
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it would not be possible to distinguish such an event from a single pN colli- 
sion with the mass Mi + M 2 produced in the center. Therefore hadron-nucleus 
and nucleus-nucleus collisions, rigorously speaking, only give information on the 
events with two fixed rapidity gaps, without specifying the exact nature of the 
object produced in between. 

However the situation improves if one takes into account that the CD events 
have a very small probability, so that the described event with two central masses 
produced is highly improbable. Then the typical CD event will be that the central 
mass is produced in only one of the inelastic collisions, all others belonging to 
elastic or low mass diffractive events. In such a case the M^-dependence of the 
CD cross section will evidently repeat that for pp collisions. 

An appropriate tool to calculate the CD cross section in hadron-nucleus and 
nucleus collisions is the so-called "criterion C" H58, |59, |60, 61|. It has been known 



since long ago that in the Glauber model the inelastic cross section is screened 
only by itself. As it turns out, there are many other types of events which are 
screened by themselves. These events have to possess a certain property C which 
satisfies the following requirement: Any superposition of NN events, in which at 
least one has property C, leads to the hA or AB event with property C . As a 
consequence, the only possibility to obtain an hA or AB event without property 
C should be that all NN events do not possess this property. 

To translate the criterion C into formulas, let <J l ^ B {a) be the inelastic nucleus 
A-nucleus B cross section considered as a function of the total nucleon-nucleon 
cross section a. Then the criterion C tells that for the events with the property 
C 

°AB = (21) 

which exactly means that they are shadowed only by themselves. 

Passing to CD events, we have to apply the criterion C twice. First choose 
as events satisfying the criterion C those with particles produced at least in one 
of the two fixed rapidity gaps which determine the central region. Evidently any 
superposition of such NN events leads to an AB event of the same type. The 
complementary events are those in which no particle is produced in any of the 
two gaps. For NN collisions they include elastic events plus events in which 
particles are also produced in between the gaps and in the low mass diffractive 
regions of the projectile and target, above the upper gap and below the lower gap. 
Therefore the nucleon-nucleon cross section a c will be given by the difference 

a = a — a —a —a \ZZ) 

where a lmd refers to the mentioned low-mass diffractive contribution. a cd denotes 
the central diffraction contribution (to which also particles scattered elastically 
or low-mass diffraction contribute). Note that the magnitude of the low-mass 
diffractive part depends on the choice of the experimental setup which defines 
the events satisfying the criterion C. In particular, if these events are chosen to 
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include all low-mass diffractive contribution to the NN cross section, then the 
term a will not appear in (ffif) and the cross section a cd will include only events 
with some particle produced in between the gap plus the projectile and target 
nucleons scattered elastically. In the following we shall be interested precisely in 
such CD events, having in mind that in the experimental setup for the photon- 
photon interaction small values of the momentum transferred to the colliding 
nuclei dominate, for which excitation of nucleon resonances is prohibited. For 
that reason we assume that all events with excitation of nucleon resonances are 
included in the events with the property C. Then (|22"D simplifies to 

_c An ^cd / OQ\ 

a = a —a [26 ) 

Subtracting from the total AB cross section the one with the described prop- 
erty C we find a cross section which is a sum of the elastic cross section, the cross 
section for the diffractive dissociation of the colliding nuclei and the CD cross 
section in which the nuclei either stay intact or are diffractively dissociated 

,_eZ i difd | cd ^tot An I An ~.cd\ i r> A\ 

Vab + Vab + °ab = °ab- °ab\° ~° ) ( 24 ) 

Now we repeat this argument taking for events satisfying C those in which any 
particle is produced in the whole rapidity range spanned by the two gaps and the 
rapidity interval in between. Such events also satisfy the conditions implied by 
the criterion C . The complementary events are now those in which no particle is 
produced in the described rapidity range at all. For NN collisions these are pure 
elastic events. For AB they also include nuclei diffraction dissociation events. 
Similar to (^4j) we then get the well known formula 

el | difd tot in /_in\ /or-\ 

Oab + °ab =(J ab-Vab\v ) (25) 



Subtracting (|25|) from ( f24|) we obtain the desired CD cross section for AB 
collisions 

cd in /_in\ Jin (in „cd\ (oa\ 

°AB = °AB\° )-°AB\° -° ) (26) 

This general formula can be written in an explicit form for HA collisions where 
the explicit dependence cr^(cr) is known. In the Glauber model, for fixed impact 
parameter B, 

= 1 " (1 " °T{B)) A ~ 1 - exp(-AxT(5)) (27) 

where T(B) is the nuclear profile function. We then find the CD cross section 

ofa = (1 - (a in - a cd )T(B)) A - (1 - a in T(B)) A 

~ exp{-Aa in T{B)){exp{Aa cd T{B) - 1) (28) 

From the derivation it is clear that this cross section refers to the total prob- 
ability to produce some particles between the two gaps and does not specify the 
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particles energy M. However, as mentioned, we can make use of the fact that o cd 
is very small. Then one approximately finds from (|28[) 

af A ~ a cd AT(B)(l-a in T(B)) A ' 1 

~ <7^AT(fl)exp(-A<7 in T(B)) (29) 

This cross section is linear in a cd . Therefore we can easily find thep distribution 
in M 



dM dM 
da cd (M) 
dM 



AT(B) exp(—Aa in T(B)) (30) 



As we observe, the absorption factor is universal and does not depend on M. 

The application of this formalism to nucleus-nucleus collisions is hampered by 
the absence of an explicit (and tractable) expression for cr™ B (cr). We shall use the 
well-known optical approximation in which (T l AB (a) is given by the same formula 
( |27|) with A — > AB and an effective profile function for the two colliding nuclei 

Tab = J d 2 BiT A (Bi)T B (B - B x ) (31) 

where Bi is the usual two dimensional impact parameter. Our final formula for 
the CD cross section in AA collisions then follows from (j30|) 

da cd A (M) _ da cd {M) a 

dM ~ dM AeS {6Z) 

where Al s is the "effective" atomic number of the colliding nuclei defined by 

A 2 cS = A 2 J d 2 BTaa(B) exp(-A 2 a m T AA (B)) (33) 
We compute A 2 S using Woods-Saxon nuclear densities 

p( r ) = ^—rr (34) 

rV ' 14- e (r-r )/a ^ > 

with the standard parameter values [B2| r = 1.14 A 1 / 3 fm and a = 0.545 fm and 
and o m = 73 mb, as predicted in |29[ for ^Js = 6 TeV. The values of A 2 S are 
given in Table [l] . 

One observes that they are much smaller than A 2 . Actually A 2 S ~ A 1 / 3 , so 
that even for very heavy colliding nuclei the CD cross section is only an order of 
magnitude greater than for proton-proton collisions. 
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2.3.2 Single diffraction cross sections in hadron-nucleus and photon- 
nucleus collisions 



The derivation given above can be applied also to single diffraction with only 
one rapidity gap from the projectile or target side. One only has to appropri- 
ately change the events satisfying the criterion C . Then, instead of the central 
diffractive nucleon-nucleon cross section a cd in Eqs. ( p2"D - fl3"0|) and (|32|), the 
corresponding single diffractive cross sections a d will appear. In particular, the 
final result for the nucleus-nucleus single diffractive cross section will be 

d in ( „in\ in ( in /nr\ 

® AB = °ABK° )-°AB\P ~ ° ) ( 35 ) 

For hA collisions we obtain instead of (28) and (29) 



a d hA = ex V {-Aa m T{B)){ex V {Aa d T{B) - 1) 

~ a d AT(B)exp{-Aa in T{B)) (36) 

and an analogous formula for da d A /dM 2 similar to (|30|). 

All of these formulae can be applied also for jlP interactions in heavy ion 
collisions. 

The final formula for jA diffractive scattering is similar to (|32|) 

jtf) d^hiM) 

with 

A cS = Aj d 2 BT A (B) exp(-Aa m T A (B)) (38) 



We calculate A eS with the same input as for A 2 eS in Eq. Q31f ) and present the 
values also in Table |l| Actually also A c g ~ A 1/3 , this is a behaviour found before 
by Ranft and Roesler pB[ and by Faessler pi . 



3 Comparing hadronic photon-photon interac- 
tions with diffractive interactions in heavy ion 
collisions 

3.1 Cross sections 

We calculate the cross sections da/dMx for the production of a central cluster 
of hadrons with invariant mass Mx- This is done for proton-proton collisions as 
well as for symmetrical heavy ion collisions with the same projectile and target 
ion. The calculations are done for light (O), medium (Ca, Fe, Ag) and heavy 
ions (Pb) with the energies of the future LHC hadron and heavy ion collider 
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under construction at CERN. The two ions most likely to be used in the LHC 
experiments 0] are Pb for obtaining the highest energy densities and Ca, where 
a larger luminosity than with Pb can be obtained. It was shown that the 
luminosity in Ca-Ca collisions could be up to a factor 10 4 larger than in Pb-Pb 
collisions. This gain more then compensates the loss in the two-photon flux in 
the Ca-Ca reaction. (L Pb _ Pb ?s 5 x 10 26 cm~ 2 s _1 , L Ca _ Ca ?s 5 x 10 30 cm" 2 s _1 ). 

For photon-photon collisions we use two different approaches for calculating 
photon fluxes described in the Appendix. These two approximations give slightly 
different cross sections especially at large masses Mx- 

In Figs. f|-^ we compare the cross section for the production of a hadronic 
cluster of invariant mass Mx in photon-photon collisions with the correspond- 
ing cross section for the double-pomeron and photon-pomeron reactions. The 
calculation was done for pp collisions and for heavy ion collisions O-O, Ca-Ca, 
Fe-Fe, Ag-Ag, and Pb-Pb at the energy of the LHC (y/s = 6 A TeV). For the 
photon-photon collisions we show the results using the form factor approximation 
(F) and the semi-classical, geometric approximation (G) to calculate the heavy 
ion photon flux. For the Pb-Pb reaction we use in fact two different form fac- 
tors, (F) the geometric approximation to the form factor (see Eq. (|47D ) and (FF), 
the Gaussian approximation. Both approximations lead to nearly identical cross 
sections da/dMx, therefore, we present for all other reactions only the (F) cross 
section. The double-pomeron cross sections are given for three different kinemat- 
ical cuts (M cd > 2 GeV/c 2 , c = 0.90, 0.95 and 0.97). The photon-pomeron cross 
sections are given for only one setting of the kinematical cuts (M 7 jp > 2 GeV/c 2 
and c = 0.95). Of course, these cuts depend on the experimental setup of such 
an experiment, which is not known at present, but any cuts are easy to apply to 
our Monte Carlo events. 

Even given the considerable uncertainty in the cross sections for CD and 
photon-pomeron collisions, the conclusions from Figs. |]-|] are rather obvious. The 
study of the pure two-photon reaction without background from central diffrac- 
tion and photon-pomeron collisions is not possible. In proton-proton collisions 
and collisions of light and medium heavy ions the central diffraction reaction dom- 
inates. For heavy ion collisions the photon-pomeron cross section is comparable 
to the photon-photon cross section. Many interesting particles might be produced 
and studied in photon-photon photon-pomeron as well as pomeron-pomeron col- 
lisions. If a reaction is to be studied using pomeron-pomeron collisions, then the 
best results should be obtained in pp collisions, where the highest luminosity can 
be obtained. 



3.2 Rapidity distributions 

In order to demonstrate that the three reactions studied lead to very similar 
signatures of the events we present in Fig. |10| pseudorapidity distributions of 
the produced hadrons. The Figure shows the distribution of hadrons in photon- 
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photon reactions in Pb-Pb and Ca-Ca heavy ion collisions (only the hadrons 
produced in the central cluster of particles are included in the histogram). The 
distribution in CD are presented for pp collisions and in this case we include also 
the scattered incoming protons into the histogram. The distribution for photon- 
pomeron collisions is presented for Pb-Pb collisions and includes only the hadrons 
from the central cluster. 

With the photon fluxes used and with the kinematical cuts for the diffractive 
reactions, we obtain in all cases two large rapidity gaps between the central cluster 
of hadrons and the scattered original protons or heavy ions. The distributions 
given in Fig. [1C] represent the average pseudorapidity distributions averaged over 
many collisions and the mass spectrum given in Figs. (| to || In each single 
event the pseudorapidity distribution corresponding to the central cluster will 
be less wide than the average and it will in general not be in the center of the 
nucleon-nucleon CMS. 

4 Conclusions and summary 

Photon-photon collisions as well as central diffraction and photon-pomeron col- 
lisions are very interesting reaction channels at proton-proton and heavy ion 
colliders in the TeV energy range. 

Here we compare the cross sections of the three channels for the production of 
hadronic systems of given invariant mass Mx- The three reaction channels have 
very similar experimental signatures, a central cluster of produced particles and 
two large rapidity gaps. 

The cross section for CD and photon-pomeron collisions is still subject to 
large uncertainties (as large as a factor of three) in the TeV energy region. 

We find for proton-proton and light up to medium-heavy ion reactions that the 
central diffraction cross section dominates the two photon and photon-pomeron 
cross sections. In collisions of the heaviest ions like Pb-Pb the photon-pomeron 
channel and the photon-photon channel are of comparable magnitude, however 
at very large masses CD can be larger than the other two channels depending on 
the experimental cuts. 

The conclusion is obvious. Central double-pomeron processes are best studied 
in pp collisions, where the largest luminosities are obtained, this compensates the 
rise of the double-pomeron cross section in heavy ion reactions. If photon-photon 
reactions are to be studied, then it should be the best to use the heaviest ions 
available at future colliders but even then the background from photon-pomeron 
collisions will be roughly of the same size as the signal. 
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A Appendix:Photon flux calculation 

It is convenient to define the luminosity function for the photon flux in hadron- 
hadron scattering by 

f{yi,Pl;V2,pl)®{si,2 - Snrin). (39) 



dy x dP{dy 2 dPi 

Here, p\ and p 2 are the four momenta of the photons forming the subsystem. 
p\ = —P 2 and p\ = — P| are the photon virtualities. The variables y\ and y 2 
denote approximately the energy fractions taken by the photons from the initial 
hadrons as explained below. The Heavyside function in ( |39"1) restricts the invariant 
mass of the system formed by the four momenta p± and p 2 to allow the application 
of the model. To obtain the cross sections for the processes mentioned above, the 
luminosity function is folded with the 77 cross section. 

As a first step, photon emission off a pointlike particle is discussed using 
as example the kinematics of ep scattering, shown in Fig. [TT|. To characterize 
deep-inelastic scattering, we use the variables x and y 

x = ^ y = {P ' Pp) (40) 

2 (Pp • P) (p e ■ Pp) 

where x denotes Bjorken's scaling variable. Then, the differential cross section 
for ep scattering via photon exchange can be written in terms of the structure 
functions Fi(x, P 2 ) and F 2 (x,P 2 ) 

4:ira 2 m \ / 2m 2 \ ^ , ~ 

J zi/(i--prJ F i& p ) 

p « F 2 x,P 2 ) . 




((Pe-Pp) 2 - ml - ml) 2 



(41) 

Here, a em denotes the fine structure constant and m e is the electron mass. In the 
limit of high collision energies, terms proportional to m 2 / (p e -p P ) 2 and m 2 / (p e -p P ) 2 
can be neglected. Using the optical theorem, the structure functions can be re- 
lated to the total jp cross sections for virtual photons with transverse polarization 



16 



(helicity ±1) <jq? and scalar polarization (helicity 0) <jg P |66 



F 1 (x,P 2 ) = fa^K 1 *) a v (42) 
F(xP 2 ) = 2 P 2 (l-*){°? + °J) m 

For small values of P 2 , simplifies to 

- -2=L ( 1 + (1 ~^ )2 - 2m 2 y — ^ aZ p (44) 

It should be emphasized that in the high-energy limit the flux of weakly virtual 
photons can be factorized out and is independent of the second scattering particle, 
which allows to introduce the generic photon flux function for bremsstrahlung 



P 2 ) = 7^rA " l ~ 2mly-^ . (45) 




Due to the complex structure of the charge distribution p{x) in hadrons, 
several approximations and assumptions are necessary to calculate the flux of 



weakly virtual photons. The approaches in literature can be subdivided into [g7 
(i) methods using charge form factors for the hadrons |2], |5| and (ii) methods 
using the semi-classical methods and geometrical interpretations on the basis of 
the impact parameter representation (for example || p], |6^, |j ) . 

In the form factor approach, Eq. fl45|) can be used directly. The effects due to 
the finite charge space-distribution can be included by substituting 

«em-^ Z 2 a cm \F(p 2 )\ 2 (46) 

for each colliding hadron where Z denotes the electric charge number. The weak 
point on this approach are the almost unknown elastic form factors F(p 2 ) for 
heavy ions. The simplest assumption for the heavy ion elastic form factor is 
motivated by the geometrical interpretation: In the classical picture one should 
only consider photons having an impact parameter B relative to the hadron 
greater than the transverse hadron size R ~ 1.2fm A 1 / 3 . With P 2 ~ 1/B 2 
follows 

f 1 , -p 2 < 1/R 2 
F(p 2 ) = / d 3 x p(x) e^ s = I (47) 

1 , -p 2 > 1/R 2 

More realistic parametrizations of the elastic form factor can be found in literature 
3S]. For example, the form factor of 206 Pb can be parametrized by a Gaussian 
distribution F(p 2 ) = exp(p 2 /Ql) with Q « 55 — 60 MeV as used in Ref. ||. 

The basis of the semi-classical (geometrical) methods is the fact that a fast- 
moving charged particle develops a magnetic field almost of the same size as the 
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electric field. This can be described by photons moving parallel to the particle 
at an impact parameter B (see Fig |T2|) . 

The number of equivalent photons is given by [70] 



f(y,B) 



Z CX-em i \ 2 1 
5— \ m V) - 

7T y 



2 

m 



K 2 1 (m\B\y) + —K 2 (m\B\y) 



(48) 



where i^o and K\ denote the modified Bessel functions, E and m are the energy 
and the mass of the hadron (heavy ion), respectively. Since the virtualities of 
the photons are neglected, the photon energy is given directly by u = yE. For 
elastic heavy ion scattering, the impact parameter of the equivalent photons is 
restricted to \B\ > R. The total photon flux follows from 

f(y) = 2tt / f(y,B)BdB (49) 

JR 

-2 



2Z 2 a e 



y 7T 



(50) 



with £ = mRy. The transverse distance between the particles should be larger 
than the transverse sizes of the particles to make sure that the particles do not 
interact hadronically which leads to the luminosity function j3], |68], § 



dL 



L L f(yi,Bi)f(y 2 ,B 2 ) 

J|Bi|>i?i J\B 2 \>R 2 



dy\dy 2 i|Bi|>fli 

x0(s 1>2 - s min ) e(|J3i - B 2 \ - {Ri + R 2 )) d 2 B 1 d 2 B 2 . 



(51) 



In the case of photon-pomeron scattering one has to consider photon-hadron 
(photon-heavy ion) scattering. The flux function reads (for example, for the 
electromagnetic interaction of particle B) 

dL f 

— = /(2/2,S 2 )e(s lj2 - Smin ) . (52) 

dy 2 J\b 2 \>r 

Here, the radius R is the minimum impact parameter of the photon that the 
heavy ions do not overlap in transverse space as shown in Fig. |T3|. A value of 
R = R\ + R 2 would satisfy this condition, however, since diffractive processes are 
mainly peripheral processes p3|, there may be also contributions with R 2 < \B\ < 



Ri + R 2 . The results shown here were obtained using R = R\ + R 2 . Lowering 
the impact parameter cutoff to R = R 2 increases the photon flux approximately 
by 20%. 
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Table 1: Results for the effective A 2 e?i in AA collisions (second column) and the 
effective A-dependence in hA collisions (third column) using a in = 73 mb and 
Woods-Saxon nuclear densities. 



Nucleus 


A 




^cff 





16 


4.94 


2.36 


Ca 


40 


6.21 


2.92 


Fe 


56 


6.76 


3.16 


Ag 


108 


8.00 


3.71 


Pb 


208 


9.52 


4.39 




a) b) 



Figure 1: High- mass photon diffraction dissociation: diffractive cut of the triple- 
pomeron graph a) and the corresponding chain system b) 
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Figure 2: Central diffraction: diffractive cut of the double-pomeron graph a) 
and the corresponding chain system b). 
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Figure 3: The energy dependence of the central diffraction cross section. We 
compare the cross section as obtained from Phojet with unitarization using 
a supercritical pomeron with the cross section obtained by Streng [OT without 
unitarization and with a critical pomeron. Both cross sections are for the same 
three kinematical cuts: M cd >2GeV/c 2 and c =0.90, 0.95 and 0.97. Please note, 
to identify the cross sections, the cross sections decrease with rising c. 
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Figure 4: We compare the cross section for the production of a hadronic cluster 
of invariant mass M x via photon-photon interaction in proton-proton collisions 
using two different methods (F) and (G) to calculate the photon flux, described 
in the Appendix with the corresponding cross section for the diffractive reactions. 
The central diffraction cross sections (pomeron-pomeron collisions) are given for 
three different kinematical cuts M cd > 2 GeV/c 2 , c = 0.90, 0.95 and 0.97. The 
single diffraction photon-pomeron cross section is given for M lF > 2 GeV/c 2 and 
c = 0.95. 
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Figure 5: As Fig. [| but for heavy ion reactions O-O. 
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Figure 6: as Fig. f|but for heavy ion reactions Ca-Ca. 
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Figure 7: As Fig. |] but for heavy ion reactions Fe-Fe. 
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Figure 8: As Fig. [| but for heavy ion reactions Ag-Ag. 
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Figure 9: As Fig. [| but for heavy ion reactions Pb-Pb. 
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Figure 10: The pseudorapidity distribution of charged hadrons in photon-photon 
collisions in Pb-Pb and Ca-Ca heavy ion collisions (only the hadrons produced 
in the central cluster of particles are included in the histogram) compared to the 
corresponding distribution in central diffraction in pp collisions (in this case also 
the scattered original protons are included in the histogram) and to the corre- 
sponding distribution in photon-pomeron collisions in Pb-Pb collisions (in this 
case again only the hadrons in the central cluster are included in the histogram). 
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Figure 11: Diagram of ep scattering via one-photon exchange. 
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Figure 12: Semi-classical model of the photon-photon scattering in hadron- 
hadron interactions. 
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Figure 13: Semi-classical model of the photon-hadron scattering in hadron- 
hadron interactions. To obtain a model for photon-pomeron scattering, only 
single diffractive events where the photon dissociates are considered. 
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